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ABSTRACT 

Context Secondary positrons are produced by spallation of cosmic rays within the interstellar gas. Measurements have 
been typically expressed in terms of the positron fraction, which exhibits an increase above 10 GeV. Many scenarios 
have been proposed to explain this feature, among them some additional primary positrons originating from dark matter 
(*~^ annihilation in the Galaxy. 

Aims. The PAMELA satellite has provided high quality data that has enabled high accuracy statistical analyses to be 
made, showing that the increase in the positron fraction extends up to about 100 GeV. It is therefore of paramount 
importance to constrain theoretically the expected secondary positron flux to interpret the observations in an accurate 
way. 

^ Methods. We focus on calculating the secondary positron flux by using and comparing different up-to-date nuclear 

cross-sections and by considering an independent model of cosmic ray propagation. We carefully study the origins of 
the theoretical uncertainties in the positron flux. 

Results. We find the secondary positron flux to be reproduced well by the available observations, and to have theoretical 

uncertainties that we quantify to be as large as about one order of magnitude. We also discuss the positron fraction 
I— I issue and find that our predictions may be consistent with the data taken before PAMELA. For PAMELA data, we find 

^ <~i that an excess is probably present after considering uncertainties in the positron fiux, although its amplitude depends 

^ I strongly on the assumptions made in relation to the electron fiux. By fitting the current electron data, we show that 

when considering a soft electron spectrum, the amplitude of the excess might be far lower than usually claimed. 
Conclusions. We provide fresh insights that may help to explain the positron data with or without new physical model 
ingredients. PAMELA observations and the forthcoming AMS-02 mission will allow stronger constraints to be aplaced 
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52 on the cosmic-ray transport parameters, and are likely to reduce drastically the theoretical uncertainties. 
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00 A&A 501, 821-833 (2009). annihilation of dark matter particles. The positron excess 

Preprint LAPTH-1272/08 and DFTT 25/2008 expected in this framework is very uncertain, because the 

nature of dark matter is unknown, and the propagation 

_ of positrons involves physical quantities that currently arc 

1. Introduction qI^q jjq^ precisely known. The related astrophysical unccr- 

O Among the different particles observed in cosmic rays, ^i^^^if calculated and quantified in Delahaye et al. 

^ positrons still raise unanswered questions. Cosmic positrons (2008), where it was shown that they may be significant, 

are created by spallation reactions of cosmic ray nuclei with especially m the low energy part of the spectrum, a prop- 

^ interstellar matter and propagate in a diffusive mode, be- ^'^^ common also to the antiproton (e^g- Donato et al. 2004) 

cause of their interaction with the turbulent component of ^nd the antideuteron (Donato et al. 2000, 2008) signals. For 

, . the Galactic magnetic field. The expected flux of positrons Positions, sizeable fluxes from dark matter annihilation axe 

^ can be calculated from the observed cosmic ray nuclei tyP^caUy possible if dark matter overdensities are locally 

fluxes, using the relevant nuclear physics and solving the f ^^^^t' Y^'f" f ^^^Y into the so-called boost 

diffusion equation factor . A detailed analysis of the admissible boost fac- 

The HEAT experiment (Barwick et al. 1997; Beatty tors for positrons and antiprotons was performed by Lavalle 

et al. 2004) showed that the positron fraction (the ra- ^l' (2008b), who showed that boost factors are typically 

tio of the positron o the total electron-positron fluxes) confined to less than about a factor of 10-20. Computing 

possibly exhibits an unexpected bump in the 10 GeV re- the antimatter fluxes directly m the frame of a cosmological 

gion of the spectrum. Although this bump could be due N-body simulation leads to the same conclusions (Lavalle 

to some unknown systematic effect, the HEAT result has ®* 2008a). 

triggered many explanations. For instance. Moskalenko & The PAMELA experiment (Picozza et al. 2007) has re- 
Strong (1998) suggested that an interstellar nucleon spec- leased its first results on the positron fraction for ener- 
trum harder than that expected could explain the excess, gies ranging from 1.5 GeV to 100 GeV and with a large 
Many works also focused on the dark matter hypothesis, statistics (Adriani et al. 2008). The positron fraction is 
the bump being due to a primary contribution from the observed to rise steadily for energies above 10 GeV, re- 
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inforcing the possibihty that an excess is actually present. 
It is therefore timely and crucial to complete a novel anal- 
ysis of the positron flux, including a robust estimation of 
the accuracy of the theoretical determination. The calcula- 
tion of the uncertainties affecting the standard spallation- 
induced positron population is especially important when 
unexpected distortions are observed in experimental data, 
to address the issue in a more robust way. This and an 
analysis of the positron signal from dark matter annihi- 
lation and its astrophysical uncertainties (Delahaye et al. 
2008), will set the proper basis for discussing in detail the 
experimental results. 

The uncertainties in the positron flux have several ori- 
gins. First, the cosmic ray nuclei measurements have their 
experimental uncertainties, which then affect the predic- 
tions of induced secondary fluxes, such as positrons. Second, 
various modelings of the nuclear cross-sections involved in 
the positron production mechanism are available, and they 
are not in complete agreement with each other, implying 
a range of theoretical variation. Third, the uncertainties in 
the propagation parameters involved in the diffusion equa- 
tion were thoroughly studied (Maurin et al. 2001). A de- 
tailed analysis of their impact on the secondary positron 
flux is therefore needed. 

In this paper, we therefore study the secondary positron 
production and transport in the Galaxy, with emphasis on 
determining the various sources of uncertainties, namely: (i) 
the nuclear cross^ections, (ii) the effect induced by the pri- 
mary injection spectra, (iii) the local interstellar medium, 
and (iv) the propagation modeling. The paper is organized 
into three main parts. The positron injection spectrum and 
its uncertainties are derived in Sect. 2. We use up-to-date 
nuclear cross-sections and show the differences with older 
parameterizations. The Green functions associated with the 
positron propagation throughout the Milky Way are dis- 
cussed in Sect. 3. Our slab model is mostly characterized 
by energy losses and diffusion caused by magnetic turbu- 
lences. The secondary positron flux at the Earth is pre- 
sented in Sect. 4 with a range of variation that includes 
the effects discussed throughout the paper. We also con- 
firm that considering the cosmic ray proton and a retro- 
propagation as well as diffusive reacceleration and convec- 
tion has little effect on our results above a few GeV. Our 
results for the positron flux and its uncertainties agree with 
all the available measurements by different experimental 
collaborations (we recall that PAMELA does not provide, 
at the moment, the positron flux, but only the positron 
fraction). The positron flux, by itself, does not exhibit un- 
usual features and is in basic agreement with the results of 
Moskalenko & Strong (1998) and Porter et al. (2008). This 
could imply that both the HEAT excess and the PAMELA 
rise observed in the positron fraction, originate from dark 
matter annihilation, but we argue in Sect. 5 that electrons 
might also play an important role. Depending on the data 
set used to constrain the electron spectrum, the estimate of 
the positron fraction can exhibit quite different behaviors, 
which can affect the interpretation of the positron fraction 
data. We show that our predictions can indeed be consistent 
with the existing data of the positron fraction for a soft elec- 
tron spectrum, compatible with flts to experimental data. 
Hard electron spectra, instead, definitely point toward the 
presence of an excess. We finally insist in Sect. 6, as a main 
result of our analysis, that the current uncertainties in the 



transport parameters translate into one order of magnitude 
uncertainties in the secondary positron spectrum. 

2. Production of positrons by spallation 

Secondary positrons are created by spallation of cosmic 
ray nuclei (mainly protons and helium nuclei) on inter- 
stellar matter (mainly hydrogen and helium). We compute 
qe+(x,_Be), the number of positrons of energy E^^ created 
per unit volume at position x, per unit time and per GeV. 
The positron source term reads: 

ge+(x,£;e) = 47r ^ ^ ntarg(x) (1) 

targ— H,Hc proj— p.o: 

/(Jb(7 
^proj (x, £'proj) X dEp^fjj X -^-^{Eproi ~^ -^e)) 

where $proj (x,-Eproj) denotes the cosmic ray nucleon flux 
at position x, ntarg(x) the number density of target nu- 
clei, and da/dEg the cross-section for the reactions creating 
positrons. We discuss these quantities this section. 

2.1. Spallation cross-sections forp + p-^e^ 

We flrst focus on the differential cross-section for the pro- 
duction of positrons. This production occurs by means 
of a nuclear reaction between two colliding nuclei, yield- 
ing mainly charged pions tt^ and other mesons, for which 
positrons are one of the flnal products of the decay chain. 
There are four main possible collisions: cosmic ray (CR) 
proton on interstellar (IS) hydrogen or helium; CR alpha 
particle on IS proton or helium. For the sake of clarity, we 
present only the formulae for the proton-proton collisions, 
but include all four processes in our results on the positron 
spectra. 

At energies below about 3 GeV, the main channel for 
production of positrons involves the excitation of a Delta 
resonance, which then decays into pions: 

p + H-p + A+-{P + ''°+ . (2) 
[n -I- 7r+ 

The charged pions decay into muons, which subsequently 
decay into positrons. At higher energies direct production 
of charged pions proceeds with the process: 

p-hH^p-Fn-h7r+ . (3) 

Kaons may also be produced: 

p + H^X + ii:± , (4) 

and the decay of kaons produces muons (63.44 %) and pions 
(20.92 %), which then decay into positrons as final products 
of their decay chain. 

To compute the differential cross-section for the pion- 

production processes, we need the probability da{Ep 
E^ ) I dEy^ of a spallation of a proton of energy Ep yielding 
a pion with energy E^^ and the probability V^E-n Eg) 
of such a pion eventually decaying into a positron of en- 
ergy Eg. This second quantity can be computed thanks to 
basic quantum electrodynamics, whereas several parame- 
terizations of the first quantity can be found in Badhwar & 
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Fig. 1. Comparison between various parameterizations of 
proton energies. 
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the positron production cross-section at different incident 



Stephens (1977), Tan & Ng (1983) and Kamae et al. (2006). 
The production cross-section of positrons is then given by: 

da , , 
— (i^p^.-^i..)^ 

■^{Ep ^ E^) X dE^ X V{E^ ^ Ee) . (5) 

Kamae et al. (2006) also provided a direct parameteriza- 
tion of the p + p ^ reaction ^ . All the afore mentioned 
parameterizations differ from each other and have been cal- 
ibrated with different nuclear data sets. The default choice 
in our calculations, unless stated otherwise, is the Kamae 
et al. (2006) parameterization, which includes additional 
processes (especially resonances other than the Delta at low 
interaction energies) and has been calibrated with recent 
data. Ss stated in the original paper, we nevertheless warn 
that this parameterization relies on fits to Monte Carlo sim- 
ulations and may be easilly affected by small uncertainties. 

The cross-section for the process involving kaons can 
be computed in a similar way to the calculation of direct 
pion production. The QED expressions for the production 
of positrons from the kaon are provided e.g. in Appendix D 
of Moskalenko & Strong (1998). 

In Fig. 1, we plot the cross-section for the positron pro- 
duction from the p-H scattering, as a function of positron 
energy. The incident proton energy is set to 2, 10 (left), 
50 and 200 (right) GeV. The three different plots at fixed 
proton energy correspond to the cross-section parameteri- 
zations of Kamae et al. (2006) (solid), Badhwar & Stephens 
(1977) (dashed) and Tan & Ng (1983) (dotted). The differ- 
ences between these plots vary with both incident proton 



^ A few typos in the published version have been corrected 
with the kind help of the authors. 



and final positron energies. For protons of intermediate en- 
ergies of 10 and 50 GeV, the flatter part of the cross-section 
(at low energies) varies only slightly between the different 
parameterizations. The Kamae et al. (2006) parameteriza- 
tion undershoots and then overshoots the other two models 
only in the high energy range, while for slow protons (see 
the 2 GeV case), it provides more positrons. Because both 
the multiple baryonic resonances around 1600 MeV and the 
standard A(1232) state have been added into their model. 

2.2. Incident proton flux 

The proton flux <i>p (x, Ep) has been measured at the loca- 
tion of the Earth x = X0 by several experiments. Various 
parameterizations of these measurements are found in the 
literature. In this analysis, we adopt the determinations of 
Shikaze et al. (2007) and Donato et al. (2001). The effect 
induced on the positron source term in Eq. (1) is displayed 
in Fig. 2, where we also show the effect arising from the 
different nuclear physics parameterizations discussed in the 
previous subsection. The solid lines refer to Kamae et al. 
(2006), the dashed lines to Badhwar & Stephens (1977), 
and the dotted lines to Tan & Ng (1983). For each set of 
curves, the thick lines are obtained for the proton spec- 
trum parameterization of Donato et al. (2001), while the 
thin lines refer to Shikaze et al. (2007). Figure 2 illustrates 
that different parameterizations of the incident proton flux 
cause less signiflcant uncertainty than the nuclear models. 
At low (1 GeV) and high (10'^ GeV) energies, the Shikaze 
et al. (2007) and Donato et al. (2001) fluxes produce almost 
identical results on the proton flux, but at intermediate en- 
ergies the differences are also neghgible. 

In contrast, the cross-section parameterizations pro- 
duce more significant variations. Figure 2 shows that the 
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Fig. 2. Comparison of the effect due to different parameter- 
izations for the cosmic ray proton spectra on the positron 
source term, as a function of the positrons energy. The ad- 
ditional effect induced by the different nuclear physics pa- 
rameterizations is also shown. The galactic protons density 
is taken at 1 hydrogen atom per cm^. 



model of Kamae et al. (2006) gives the faintest spectrum 
in all the energy range (from 1 GeV to 1 TeV), while that 
of Tan & Ng (1983) corresponds to the maximal spectrum. 
The most significant difference occurs at 20-30 GeV, of 
about a factor of two. 

The spatial dependence of the proton flux is determined 

by solving the diffusion equation and normalizing the re- 
sulting flux to the solar value. This procedure will be re- 
ferred to as retro-propagation in the following. We find that 
retropropagation is not crucial and that one can safely ap- 
proximate the proton fiux as being homogeneous and equal 
to the solar value. This is because most positrons detected 
in the solar neighborhood have been created locally, over a 
region where the proton fiux does not vary significally. This 
will be discussed fTirther in the final section. 

For the sake of simplicity, we show only the impact of 
the pH interaction on the determination of the positron 
spectrum, but we recall that our results in the subsequent 
sections refer to cosmic ray protons and a interacting with 
both IS hydrogen and helium with densities of nn = 0.9 
cm~^ and riHe = 0.1 cm~^ respectively. These average val- 
ues for the hydrogen and helium densities are of course 
approximations based on local estimates, but we do not 
expect significant changes in the averaging within the kpc 
scale (Ferricre et al. 2007), whic;h, as we show in the sub- 
sequent sections, is the relevant scale for the secondary 
positron propagation. Cross-sections for heavy nuclei were 
dealt with in the same way as in Norbury & Townsend 
(2007). 



3. Propagation 

In the Galaxy, a charged particle travelling between its 
source and the solar neighborhood is affected by several 
processes. Scattering by magnetic fields leads to a random 
walk in both real space (diffusion) and momentum space 
(diffusive reacceleration) . Particles may also be spatially 
convected away by the galactic wind (which induces 
adiabatic losses), and lose energy as they interact with 
either interstellar matter or the electromagnetic field 
and radiation of the Galaxy (by synchrotron radiation 
and inverse Compton processes). Above a few GeV, the 
propagation of positrons in the Milky Way is dominated 
by space diffusion and energy losses. 

In this paper, the diffusion coefficient is assumed to 
be homogeneous and isotropic, with a dependence on en- 
ergy given by K{E) = l3Kf){TZ/l GV)"^, where the mag- 
netic rigidity TZ is related to the momentum p and elec- 
tric charge Ze hy TZ = pc/Ze. Cosmic rays are confined 
within a cylindrical diffusive halo of radius i? = 20 kpc 
and height 2L, their density vanishing at the boundaries 
iV(|z| = L,r) = N{z,r = R) = 0. As discussed below, the 
radial boundary has a negligible effect on the density of 
positrons in the Solar System. In this section, we do not 
consider the possibility of Galactic convection and diffu- 
sive reacceleration. These processes are taken into account 
in Sect. 4.2, where we demostrate that they have little ef- 
fect. The convective wind is assumed to carry cosmic rays 
away from the Milky Way disk in the z direction at a con- 
stant velocity Vc- Diffusive reacceleration depends on the 
velocity Va of the Alfven waves. The free model param- 
eters are therefore the size L of the diffusive halo, both 
the normalization Kq and spectral index 6 of the diffu- 
sion coefficient, the convective wind velocity Vc, and the 
Alfven velocity Va (see Sect. 3.3 for additional details). 
This model has been consistently used in several studies to 
constrain the propagation parameters (Maurin et al. 2001, 
2002; Donato et al. 2002) and examine their consequences 
(Taillet & Maurin 2003; Maurin fc Taillet 2003) for the 
standard p fiux (Donato et al. 2001), the exotic p and d 
fiuxes (Maurin et al. 2004, 2006; Donato et al. 2004; Barrau 
et al. 2002, 2005; Bringmann & Salati 2007), and also for 
positrons (Lavalle et al. 2008b; Delahaye et al. 2008). The 
reader is referred to Maurin et al. (2001) for a more detailed 
presentation and motivation of the framework. 



3.1. The Green function for positrons 

The propagation of positrons differs from that of the nu- 
clei in several respects. Although space diffusion is an 
essential ingredient common to all cosmic ray species, 
positrons suffer mostly from inverse Compton and syn- 
chrotron energy losses, e.g. Moskalenko &: Strong 1998, 
whereas (anti )protons are mostly sensitive to the galactic 
wind and the nuclear interactions as they cross the Milky 
Way disk. As a result, a positron line at source leads to 
an extended spectrum after its propagation. This disagrees 
with studies of most nuclear species for which, as first ap- 
proximation, energy losses can be neglected. Consequently, 
the diffusion equation that leads to the positron number 
density N per unit of volume and energy, with the source 
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term qf,+ {x,E), becomes: 

s 



E 

Eq 



d UE 
AN + — < — N 
dE\dt 



(7e+(x,£;). (6) 



The first term is simply the diffusion coefficient written 
as Kin) « Ko{E/Eoy, where £^0 = 1 GeV is used to 
maintain the correct units throughout the paper. 
The synchrotron and inverse Compton lossc^s can be written 
as dE/dt{E) = —E'^/{EqTe)- Defining a pseudo-time: 



{E/E^) 

t = Te 



(5-1 



1-5 

and applying the following rescaling: 



(7) 



N = {E/EofN and q,+ ix,E) = {E/Eof^'qe+i^,E), 



the diffusion equation can be rewritten as: 

ON 



di 



(8) 



which is formally identical to the well-known time- 
dependent diffusion equation (Bulanov & Dogel 1974; 
Baltz & Edsjo 1999). 

It proves convenient to separate diffusion along the ra- 
dial and vertical directions. Considering a source located 
at {x,y, z,tE^) and detected at (i?©, 0, 0, t^j^), the corre- 
sponding flux depends only on the radial relative distance 
r = \rs — ro\, the distance of the source from the plane 
z = zs and the relative pseudo-time t = tEo —tss- When 
the radial boundary is taken into account, one has to use 
an expansion over Bessel functions (Bulanov & Dogel 1974; 
Delahaye et al. 2008). However, in most situations it is 
safe to ignore the radial boundary. The Green function 
Qq (r, z, t) of Eq. (8) is then given by: 



^0(-r,r-, z) 



>(f) 



exp 



4ifnf 



xg'^{z,f). (9) 



The radial behavior of the positron Green function encour- 
ages us to define the characteristic diffusion length: 



Ad = \/4ifo7 



(10) 



This length defines the scale of the positron sphere, i.e, the 
region where most of the positrons detected at the Earth 
are produced. It depends on the injected Es and detected 
Eo positron energies in terms of the pseudo-time difference 
f. For GeV energies, < 5 kpc, typically, which defines 
the local character of the positron origin. 

The effect of boundaries along z = ±L appears in only 
g^^{z,T). For convergence properties, two distinct regimes 
are worth considering (Lavalle et al. 2007) 

1. When the extension Ad of the positron sphere is smaller 
than the half thickness L of the diffusive halo, it is most 
appropriate to use the so-called electrical image formula 
(e.g. Baltz & Edsjo 1999): 



where Zn = 2Ln + (—1)" z; 



(11) 



2. In the opposite situation, a more suitable expres- 
sion is based on an analogy with the solution to the 
Schrodinger equation in an infinitely deep square poten- 
tial: expansion of the solution over the eigenfunctions of 
the Laplacian operator: 



1 



n=l 



(12) 



where: 



1\ TT 



(j)n{z) = sin[fc„(i - |^;|)] ; kn = yn - ^ (even) 

TT 

(l)'^{z) = sin \k'^{L - z)] and k'^ = n- (odd). 

The true positron propagator for a monochromatic 
point source is related to Eq. (9) by 

(E ^Es,r,z) = ^ X ^Q(f = tE - tEs ,r,z). (13) 

The secondary positron flux at the Earth is then given by 
(considering the Earth as the origin of the coordinate sys- 
tem): 



b%{E)=^^^.j^dEs. 



(14) 



X / rf^xs g^'' {E^Es,rs, Zg) qe+ (xs, Es) , 

J slab 



where \xs\ 



Zg is the squared distance from the 



source to the Earth, and /3 is the positron velocity in units 
of the velocity of light. 



3.2. Analytical solution for a homogeneous source term 

The secondary positrons originate from spallation processes 
of cosmic rays off the interstellar gas, the latter being lo- 
cated mainly inside the galactic disk. If we assume that 
the cosmic ray fluxes and the gas are homogeneous inside 
the disk and place the radial boundaries at infinity (see 
Sect. 3.1), the spatial integral in Eq. (14) can be derived 
analytically, e.g. by implementing an infinite disk of thick- 
ness [— -Jmax, 2^max]- The flux at the Earth simplifies into 



where 



/3c 



/ dEs qe+{Es) X 
Je 



teEq 
E^ 



xr?(AD), (15) 



d^xsgQ{T = iE-iEs,r,z) . (16) 



This quantity depends only on the characteristic scale Ad 
and can be expressed as 



^(Ad) 



n=—oo 



Erf 



-Erf 



„+l COS (fcn(L- ^max)) ^ ^-fe^AnV* 



(17) 
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The upper line corresponds to the electrical image solu- 
tion, while the lower line is obtained by expanding the 
solution over the eigenfunctions of the Laplacian opera- 
tor (see Sect. 3.1). Regarding the former case, we define 
= ^riL + (-l)"z„,ax and zf'^ = 2nL - (-l)"z„,ax 
where Zmax is defined to be 100 pc. The energy depen- 
dence of the solutions is hidden in the propagation length 
Ad- This simplification is helpful because it can provide 
an efficient way of checking the numerical spatial integral. 
Moreover, if the source term is shown to be almost homoge- 
neous in the thin galactic disk, then this solution fully ap- 
plies. Homogeneity is certainly justified in the energy range 
considered here but could break down at TeV energies. 

In Fig. 3, ?7 is plotted as a function of the propagation 
length Ad- This parameter can be interpreted as the ratio 
of received to produced positrons. In the regime where Ad 
is much smaller than the half-thickness L of the diffusive 
halo, positron propagation occurs as if the diffusive halo 
were infinite. In this 3D limit, the propagator Q(7){t = Ie — 
Iest^tz) describes the probability of a positron detected 
at Earth originating in the location r and z and as such 
is normalized to unity. Secondary positrons are produced 
solely in the disk and 77 equals by definition the fraction of 
the positron sphere filled by the disk. We therefore expect 
this fraction to be close to 1 when Ad is smaller than z^ax as 
the positron sphere becomes small enough to be embedded 
entirely inside the Galactic disk. In the converse situation, 
77 decreases approximately as 

ry(AD » ;^n.ax) ^ A . (18) 

Both behaviors are illustrated in Fig. 3. 

3.3. Propagation parameters 

The propagation parameters 5, Kq, L, Vc, and Va are not 
measured directly. However, it is possible to determine the 
parameter sets that are consistent with the observed prop- 
erties of nuclei cosmic ray spectra, by comparing the ob- 
served spectra with the predictions of the diffusion model. 
In Maurin et al. (2001, 2002), the secondary-to-primary ra- 
tio B/C is used to place constraints on the parameter space, 
and in Donato et al. (2002) the information provided by the 
cosmic ray radioactive species is studied. For the discussion 
below, it is convenient to isolate three sets of parameters 
labeled MIN, MED, and MAX, defined in Table 1. These 
configurations are named according to the primary antipro- 
ton signal yielded by dark matter species annihilating in the 
Milky Way halo (Barrau et al. 2002; Donato et al. 2004) and 
filling completely the diffusive halo. The half-thickness L 
increases from 1 to 15 kpc between models MIN and MAX. 
Notice that in the B/C analysis of Maurin et al. (2001), L 
is shown to be correlated with the normalization Kq. Thick 
diffusive halos are associated with high diffusion coefficients 
and hence to high values of Ad- Inspired by Fig. 3, we an- 
ticipate that the secondary positron flux should have its 
lowest value in the MAX configuration and its highest for 
the MIN model. 

3.4. Energy losses 

A detailed analysis including all energy-loss processes is 
presented in Sect. 4.2. Here, we focus on the main processes 
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Fig. 3. The integral rj is plotted as a function of the ratio 
AD/zmax- Because rj can be interpreted as the the fraction 
of the positron sphere intersected by the Galactic disk, we 
infer that it should be unity for Ad z^ax- In the converse 
regime, 77 is proportional to the ratio Zmax/AD- See text for 
further details- 



Model 


S 


Ko 


L 










[kpcVMyr] 


[kpc] 


[km/s] 


[km/s] 


MIN 


0.85 


0.0016 


1 


13.5 


22.4 


MED 


0.70 


0.0112 


4 


12 


52.9 


MAX 


0.46 


0.0765 


15 


5 


117.6 



Table 1. Typical combinations of diffusion parameters 
that are compatible with the B/C analysis (Maurin et al. 
2001). As shown in Donato et al. (2004), these propagation 
models correspond respectively to minimal, medium, and 
maximal primary antiproton fluxes. 



for our energy range of interest. At energies higher than 
about 10 GeV, the most relevant energy losses are those 
caused by synchrotron radiation and inverse Compton (IC) 
scattering of the cosmic microwave background (CMB) and 
stars photons: 

2222 
-b'°'%e) = — = — + - + ^—. (19) 

TE Tgync T* TCMB 

Each energy-loss timescale t can be calculated by means 
of both the Compton cross-section and the correspond- 
ing radiation-field energy density, as detailed in detail by 
Longair (1994). These Compton processes can indeed be 
computed in the Thomson limit. In the energy domain of 
interest here for positrons, say 10-100 GeV, the condition 
7e+i?ph < 'Tie, where 7g+ is the positron Lorentz factor 
and i?ph is the photon energy, is fulfilled for all relevant 
frequencies of the considered radiation fields. 

By studying Eqs. (7) and (10), we note that te has 
a significant effect on the positron propagation length: 
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the higher the value of te, the farther is the origin of 
detected positrons from Earth. Since the Green function 
has a Gaussian dependence on the propagation scale, the 
importance of estimating accurately the te parameter is 
reinforced, even though we expect the effect to be lower for 
secondaries than for primaries, for which the source term 
may, in some cases, exhibit a far more pronounced spatial 
dependency. In the following, we briefly discuss the differ- 
ent contributions to the energy-loss timescale in the high 
energy range. The CMB contribution is the most straight- 
forward to compute. Adopting a mean CMB temperature 
of T = 2.725 K, we readily determine tcmb = 3.77 x 10^*^s. 

For the synchrotron contribution, the local value of the 
magnetic field, which equals Tgync, remains unknown. As ex- 
plained in Beck et al. (2003), there are two different meth- 
ods of measurement. The first relies on the intensity of the 
synchrotron radiation from cosmic electrons and gives a 
value ofB-4±l/iG (Prouza & Smida (2003)). However, 
this value depends on the adopted model, particularly in 
terms of the cosmic ray electron spectrum estimation. The 
second method uses the Faraday rotation measurements of 
pulsar polarized emission, and yields B ~ 1.8±0.3 /i G 
(Han et al. (2006)). The two results are inconsistent but 
Beck et al. (2003) were able to identify further uncer- 
tainty in the second method, which, if the thermal electron 
density is anticorrelated with the magnetic field strength, 
produces a revised value of B G [1.5; 4] G, and hence 
T,y„c e [2.47 X 10i<^; 1.76 x 10^^] s. 

Finally, to evaluate the contribution of IC processes in 
the interstellar radiation field (ISRF), we rely on the study 
performed by Strong et al. (2000). They estimated the local 
value of the ISRF energy density, whose uncertainty can 
be inferred from its variation in value within the 2 kpc 
around the solar position. This provides an average local 
ISRF energy density of ?7iad ~ 2 ± 1 eV cm~'^, hence tjc G 
[6.54 X 101^1.96 X lO^S] s. 

The leading term caused by IC scattering by the ISRF 
(star and dust light) is is clearly affect by the most signif- 
icant uncertainty. By combining all the contributions, we 
derive an uncertainty range for the energy-loss timescale 
of TE e [5.17 X 10^^ 1.77 X 10^6] s. As shown in Fig. 4, the 
uncertainty caused by te does not exhibit a strong energy 
dependence. This can be understood from the expression 
of the positron Green function given in Eqs. (9) and (13), 
and from the fact that the source term of secondaries is 
close to be homogeneously distributed in the thin Galactic 
disk. From this, the secondary positron fiux roughly scales 
as y^TE, as also illustrated in Fig. 4. Since te determines 
the positron propagation scale, the related uncertainties 
should be far stronger effeact on the primary contributions 
with a far greater spatial (or time) dependency, such as for 
instance nearby astrophysical or exotic point sources. In 
the following, we have adopted as a standard the value of 
TE^lX lO^s s. 

We note that Tgync f^nd r^, are not expected to re- 
main constant throughout the entire diffusive halo (Strong 
et al. (2000)). However, as we discuss extensively in Sect. 
4.1, most of the positrons detected at Earth with energies 
greater than a few GeV, are typically of local origin. This 
is illustrated in Fig. 10, where = 1 x 10^^ s, for de- 
tected positrons of energy higher than 1 GeV: more than 
75% of the signal has originated within a distance of 2 kpc; 
the higher the detected energy, the higher this percentage. 




10' 1 10 10^ 

E [GeV] 



Fig. 4. Interstellar secondary positron fiux i?^-^$g+ as 
a function of the energy at Earth, for different values of 
the energy loss timescale te- The longer the timescale, the 
larger the flux. The scaling relation cx y/rE is also re- 
ported. 



Since we are not interested in very low energy positrons, we 
can safely neglect spatial variations in the magnetic field B 
and the ISRF energy density Uiad discussed above. 

4. The positron flux and its uncertainties 

Figure 5 displays the calculated secondary positron fiux 
modulated at solar minimum along with the most recent 
experimental data. We used a Fisk potential (p = 600 MV 
as applied in Perko (1987). The MIN, MED and MAX cases 
are illustrated by the red solid, long-dashed and short- 
dashed lines, respectively, while the yellow area denotes the 
uncertainty in the propagated fiux caused by the uncer- 
tainty in the astrophysical parameters. The nuclear cross- 
section from Kamae et al. (2006) and the Shikaze et al. 
(2007) proton and helium fiuxes were used. We adopted 
the MED prediction as our reference model. In the same 
figure, we also plot the interstellar fiux. The upper long- 
dashed curve corresponds to the MED case whereas the 
slanted band indicates the uncertainty in the Galactic 
propagation parameters. The solid line shows the IS fiux 
from Moskalenko & Strong (1998). Below - 100 GeV, the 
yellow uncertainty band is delineated by the MIN and MAX 
models. As discussed at the end of the previous section, the 
MIN (MAX) set of parameters yields the highest (lowest re- 
spectively) values for the secondary positron fiux. Since we 
considered more than about 1,600 different conflgurations 
compatible with the B/C ratio (Maurin et al. 2001), other 
propagation models become important in determining the 
extremes of the uncertainty band at higher energies. The 
maximal flux at energies above 100 GeV does not corre- 
spond to any specific set of propagation parameters over 
the whole range of energies, as already noted in Delahaye 
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Fig. 5. Secondary positron flux as a function of the 
positron energy. The blue hatched band corresponds to 
the CR propagation uncertainty in the IS prediction, 
whereas the yellow strip refers to TOA fluxes. The long- 
dashed curves feature our reference model with the Kamae 
et al. (2006) parameterization of nuclear cross-sections, the 
Shikaze et al. (2007) injection proton and helium spectra 
and the MED set of propagation parameters. The MIN, 
MED and MAX propagation parameters are displayed in 
Tab. 1. Data are taken from CAPRICE (Boezio et al. 2000), 
HEAT (Barwick et al. 1997), AMS (Aguilar et al. 2007; 
Alcaraz et al. 2000) and MASS (Grimani et al. 2002). 



et al. (2008), where the case of positrons produced by dark 
matter annihilation was studied. 

From Fig. 5, we see that the variation in the propa- 
gation parameters induces an uncertainty in the positron 
flux, which reaches about one order of magnitude over the 
entire energy range considered here. It is a factor of 6 at 
1 GeV, and smoothly decreases down to a factor of 4 or 
less for energies higher than 100 GeV. The agreement with 
experimental data is quite good at all energies within the 
uncertainty band. The Moskalenko & Strong (1998) predic- 
tion of the IS secondary positron flux as parameterized by 
Baltz & Edsjo (1999) is indicated by the black solid curve, 
and hardly differs from our reference model (long-dashed 
curve and MED propagation) above a few GeV. The HEAT 
data points are in good statistical agreement with this MED 
model. 

The effects induced by different parameterizations of 
the nuclear production cross-sections and by the variation 
in the proton injection spectrum are shown in Fig. 6. In 
the left panel, we present the TOA positron fluxes calcu- 
lated from the Tan & Ng (1983) (dotted), Badhwar et al. 
(1977) (dashed), and Kamae et al. (2006) (sohd) cross- 
section models for the MED propagation scheme and the 
Shikaze et al. (2007) proton and helium injection spectra. 
The Kamae et al. (2006) model leads systematically to the 



lowest flux. For positron energies < 1 GeV, the three cross- 
section parameterizations differ by just a few percent, while 
the differences are signiflcantly larger at higher energies. 
Figure 6 translates the uncertainties in the source term qe+ 
featured in Fig. 2. Consequently, the flux obtained at 10 
GeV with the Tan & Ng (1983) or Badhwar et al. (1977) 
parameterization is a factor of 2 or 1.6, respectively, higher 
than the reference case (Kamae et al. 2006). This trend is 
conflrmed at higher energies, although the differences be- 
tween the various models are smaller above 200 GeV. 

The uncertainties caused by the proton and helium 
spectrum parameterizations are the least relevant to this 
analysis. This is demonstrated in the right panel of Fig. 6, 
where we compare the positron flux of the reference model 
(solid lines) with the flux obtained when the Shikaze et al. 
(2007) parameterization of the incident spectra was re- 
placed by that of Donato et al. (2001). The differences are 
at most 10-15% around 10 GeV, and are negligible in the 
lower and higher energy tails. 

4.1. Spatial origin of the positrons 

At every location in the Galaxy, the positron production 
by spallation is determined by the local flux of cosmic 
ray proton and helium projectiles. Their spatial distribu- 
tion <i>(x, E) was assumed to be constant and set equal to 
the value ^q{E) measured at the Solar System location. 
However, we note that these CR primaries also diffuse in 
the Milky Way, so that their flux should exhibit a spatial 
dependence. The positron source term should vary ac- 
cordingly inside the Galactic disk. The behavior of the pro- 
ton and helium fluxes with radius r can be inferred readily 
from their measured values ^q{E) once the propagation 
parameters are selected. This so-called retro-propagation 
was implemented in the original B/C analysis by Maurin 
et al. (2001). The radial variation in the proton flux is pre- 
sented in Fig. 7 for two quite different proton energies, and 
is found to be significant. This is why we questioned the 
hypothesis of a homogeneous positron production through- 
out the disk and found nevertheless that it remains viable 
in spite of the strong radial dependence of <i>(r, E). 

This is because positrons reaching the Earth were 
mostly created locally, i.e., in a region in which the pro- 
ton flux does not differ significantly from the local value 
^q{E). We evaluated the contribution to the total signal 
from a disk of radius rgourco surrounding the Earth, which 
was modeled with the source term 

'?sourco(7', E) = qe+ (r, E) x e(rsource ^ r) , (20) 

where Q{x) is the Heaviside function and r measures the 
radial distance from the Solar System. The positron flux 
yielded by g source is 0f+ (^source, -E), whose contribution to 
the total signal (E) is plotted in Fig. 8 as a function of 
''source, for several values of the positron energy E. Most of 
the positron signal originates at short distances, especially 
at high energy. At 1 TeV, more than 80% of the positrons 
are created within 1 kpc while more than half of the 100 
MeV positrons come from less than 2 kpc. Energy losses 
are indeed quite efficient. They reduce the positron horizon 
strongly as the energy increases. This is why the CR pro- 
ton and helium fluxes are close to their solar values when 
averaged over the positron horizon scale. Taking the retro- 
propagation of projectile spectra into account has therefore 
little effect on the positron flux, as is clear in Fig. 9. 
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Fig. 6. Left: TOA and IS positron spectra for three different nuclear cross-section parameterizations: Kamae (solid), 
Badhwar (dashed) and Tan & Ng (dotted). Right: TOA and IS positron spectra for two different proton fluxes: Shikaze 
(solid) and Donato (dashed). In all cases, diffusion parameters are set to the MED case of Tab. 1. 



4.2. Diffusive reacceleration and full energy losses 



line. The diffusion equation may now be expressed as: 



Space diffusion and energy losses through inverse Compton 
scattering and synchrotron emission were the only pro- 
cesses that we had considered. We had neglected many 
other mechanisms that may also affect positron propaga- 
tion. Galactic convection can sweep cosmic rays out of the 
diffusive halo and is associated with adiabatic energy losses. 
The drift of the magnetic turbulent field with respect to the 
Galactic frame with velocity Va induces both a diffusion 
in energy space and a reacceleration of particles. This so- 
called diffusive reacceleration was discussed in the original 
analysis of Moskalenko & Strong (1998), but it does not 
appear in the fitting formula proposed by Baltz & Edsjo 
(1999). Finally, bremsstrahlung and ionization should come 
into play at low energies, below ~ 1 GeV. 

We note that once all these effects are considered, there 
is no longer an analytical solution to the diffusion equa- 
tion, at least none that we know. Synchrotron and IC en- 
ergy losses occur all over the diffusive halo, whereas the 
other energy-loss mechanisms as well as diffusive reacceler- 
ation are localized within the Galactic plane. A completely 
numerical approach is always possible of course - dealing 
for instance with a realistic distribution of gas, magnetic 
fields, and interstellar radiation fields - but this has never 
been our philosophy so far. Keeping calculations as ana- 
lytical as possible has always been our guiding principle. 
We have therefore derived an approximate solution to the 
complete diffusion equation where the IC and synchrotron 
losses have been suitably renormalized and assumed to take 
place only in the disk. Obviously this is not fully correct, 
but a close inspection of the left panel of Fig. 10 and its 
solid curve indicates these values agree closely with the ref- 
erence model featured in Fig. 5 by the long-dashed blue 



V- {-Koe^VN + Yciz)N} + 



(21) 



de 



■Aoss 



(e) N - 



Ik J 



9e+(x,e) , 



where e = E/Eq and Eq — 1 GeV as defined in Sect. 3.1. 
The convective wind varies with the vertical coordinate as: 



Vc(^) 



VcUz 

-VcVLz 



if z > 
if z < 0. 



(22) 



The energy-loss term is a combination of several contribu- 
tions: 



I loss 



e 

TE 

hV- V 



c 



p 

6he 



inverse Compton and synchrotron 

adiabatic losses 
bremsstrahlung 
19. sj' ionisation. 

(23) 

The values of the constants Ki, and Ki can be found in 
Longair (1994). Diffusive reacceleration is taken care of by 
using a coefficient 



-Kb uh e 
-Ki flu < 3 In 



E 



K,, 



2 V"^ 

9 Kr, 



(24) 



The positron density N{r, z, e) is Bessel expanded. The co- 
efficients Ni(z,e) are determined all over the diffusive halo 
except for the normalizations Ni{0,e). Each of these nor- 
malizations satisfies a diffusion equation in energy space, 
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Fig. 7. Ratio of the proton flux at radius r to the so- 
lar value from Shikaze et al. (2007). In this plot, retro- 
propagation has been taken into account, and all propaga- 
tion effects of the MED configuration (i.e., convective wind, 
spallation, and diffusion) have been included. The dot refers 
to the Solar System position in the Galaxy. 



which we solve using a Crank-Nicholson semi-implicit 
method. 

The solid line of Fig. 10 considers only space diffusion 
and energy losses by IC scattering and synchrotron emis- 
sion. When these processes are supplemented by diffusive 
reacceleration, we derive the long-dashed curve with a no- 
ticeable bump at ~ 3 GeV. Below that value, positrons 
are reaccelerated and their energy spectrum is shifted to 
higher energies. Above a few GeV, IC scattering and syn- 
chrotron emission dominate over diffusive reacceleration, 
inducing a shift in the spectrum towards lower energies. 
Positrons accumulate in an energy region where energy 
losses and diffusive reacceleration compensate each other, 
hence a visible bump which is already present in the analy- 
sis by Moskalenko & Strong (1998). The short-dashed line 
is obtained by replacing diffusive reacceleration by Galactic 
convection. The wind is active at low energies, where space 
diffusion is slow. Positrons are drifted away from the Galaxy 
and their flux at the Earth is depleted. We note that diffu- 
sive reacceleration and Galactic convection were included 
separately by Lionetto et al. (2005) in their prediction of 
the positron spectrum, with the net result of either over- 
shooting (diffusive reacceleration) or undershooting (galac- 
tic convection) the data. If we now incorporate both pro- 
cesses and add the various energy-loss mechanisms, we de- 
rive the dotted curve, which also contains a bump, although 
of far smaller amplitude. The bump cannot be distinguished 
from the solid line for energies above a few GeV. This is 
the energy region where dark matter species are expected 
to distort the secondary spectrum and a calculation based 



Fig. 8. Fraction of the positrons detected at the Earth 
which are produced within a disk of radius rgource- The 
larger the energy, the closer the source. 
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Fig. 9. Ratio of the positron flux computed with and with- 
out retro-propagation, as a function of positron energy. 



solely on space diffusion and energy losses from IC scatter- 
ing and synchrotron emission is perfectly safe. 

Below a few GeV, the situation becomes more compli- 
cated, several effects at stake modifying the blue hatched 
IS and yellow TOA uncertainty intervals in Fig. 5 as dis- 
played in the right panel of Fig. 10. Reproducing the GeV 
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Fig. 11. Electron flux parameterization. 



bump of tlie IS flux with the data now requires a higher 
Fisk potential of 850 MV. The agreement seems reasonable 
below a few GeV, although a more detailed investigation 
would require a refined solar modulation model. 



5. The positron fraction 

The question of whether a positron excess is being observed 
in cosmic ray measurements remains after many years, and 
is usually addressed in terms of the the so-called positron 
fraction, i.e., the quantity 0e+/(0e+ + (^e-)- This excess 
was pointed out by Moskalenko & Strong (1998) when they 
derived their predictions of both the secondary positron 
flux and the primary plus secondary electron fluxes. This 
led to many possible interpretations, such e.g., a poten- 
tial positron injection from dark matter annihilation in the 
Galaxy (e.g., Baltz & Edsjo 1999) or from pulsars (e.g., 
Aharonian et al. 1995; Boulares 1989; Grimani 2007). 

The PAMELA collaboration reported the positron frac- 
tion from 1.5 to 100 GeV with unprecedented statis- 
tics quality (Adriani et al. 2008). The PAMELA data 
are reported in Fig. 12, which we further explain in 
this section. Compared to the typical reference prediction 
by Moskalenko & Strong (1998), a large excess appears. Our 
new predictions of the secondary positron flux and its the- 
oretical uncertainties allow us to discuss in greater depth 
the interpretation of the excess positron fraction. We note 
flrst that we have shown in the previous sections that our 
predictions are consistent with the available positron mea- 
surements, excluding the PAMELA data since only their 
positron fraction data have so far become public. From 
these pre-PAMELA data, we therefore remark that an ex- 
cess is hardly observed when considering the positron mea- 
surements only. 

One of two crucial ingredients needed to derive the 
positron fraction is of course the electron flux. It was al- 
ready noted by Moskalenko & Strong (1998) that a change 



in the electron spectrum may affect the existence of an ex- 
cess in the positron fraction. These authors compared the 
positron fraction obtained when using their own prediction 
of the electron spectrum, with that obtained when using 
the prediction of Protheroe (1982) based on the leaky box 
propagation model, and hence illustrated the difference due 
to the electrons. Today, we can attempt to take advantage 
of the higher quality existing electron data, and comple- 
ment them with our theoretical calculation of the positron 
flux and its uncertainties. 

It is nevertheless difficult to constrain the electron flux 
using the existing data, because they are available not 
only for limited energy ranges but also exhibit some dif- 
ferences in the absolute normalizations as well as the spec- 
tral shapes. Furthermore, the measurements available are 
mainly for energies below ^ 50 GeV. For instance, the AMS 
electron data points of Aguilar et al. (2007) and Alcaraz 
et al. (2000), which are known to be among the most precise 
data sets to date, reach ~ 30 GeV only. Other complications 
are expected from astrophysical modeling arguments. For 
instance, it is likely that different spectral contributions are 
important at different energies, such as secondary electrons 
at low energy, and primary electrons from local (distant) 
astrophysical cosmic ray sources at higher (intermediate, 
respectively) energies: this may imply that a blind flt to the 
data is not useful in the absence of any modeling insights. 
Although accounting for all these subtleties is beyond the 
scope of this paper, we can illustrate the importance of the 
electron spectrum by using a simple parameterization that 
is constrained reasonably up to ~ 100 GeV, i.e., encompass- 
ing the PAMELA energy range, which fits the data at lower 
energies. We emphasize that predicting the electron flux is 
unnecessary to interpreting the positron contribution to the 
positron fraction as soon as electron data become available. 
Using the electron data directly instead of a model is also 
likely to provide a far more accurate description of any ex- 
cess. 

Making use of the electron data, we have therefore mod- 
eled an electron spectrum by a power law cx E~'^ at energies 
above a few GeV and up to 100 GeV, and allowed a scatter 
in the spectral index to account for the dispersion in the dif- 
ferent sets of data. Because the permitted range of spectral 
indices that provide good flts is signiflcant, we have applied 
restrictions by taking the central value and variance found 
by Casadei & Bindi (2004), i.e., 7 = 3.44 ± 0.03. Thus, an 
index range deflned by 7 ± 3(t allows to encompass most of 
the available data below 100 GeV. The normalization of this 
electron spectrum as well as the low energy part have been 
adapted to the AMS data, which therefore includes the so- 
lar modulation effects. The choice of AMS is motivated by 
being likely to be the experimental setup least affected by 
systematic errors, but we emphasize that other setups are 
possible. Our parameterization of the electron flux is illus- 
trated in Fig. 11, where the data from HEAT, CAPRICE, 
and AMS are presented. The dispersion observed in the 
data above a few GeV and below ~ 100 GeV, i.e., over 
the entire PAMELA energy range, is reproduced well by 
taking a spectral index of 7 = 3.44 ± 0.1, as mentioned 
above. To interpret the positron fraction measurements, 
we therefore considered two cases for the electron flux, a 
soft spectrum with index 3.53, and a hard spectrum with 
index 3.35. Finally, we note that we have also presented 
in Fig. 11 the interstellar electron model of Moskalenko & 
Strong (1998), which, while those authors have since con- 
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Fig. 10. Left panel: The reference model of Sect. 4 is featured here with various effects turned on or off. Space diffusion 
and energy losses from IC and synchrotron emission lead to the solid curve. When diffusive reacceleration is added, we get 
the long-dashed line and its spectacular bump around 3 GeV. The short-dashed curve is obtained by replacing diffusive 
reacceleration by galactic convection. The spectrum becomes depleted at low energies. Including all the processes yield 
the dotted line. Diffusive reacceleration and convection are both relevant below a few GcV and induce opposite effects. 
Right panel: The hatched blue (IS) and yellow (TOA) regions of Fig. 5 delineated by the MAX and MIN curves are 
featured here with all the effects included. Above a few GeV, we get the same results as before. Data are taken from 
CAPRICE (Boezio et al. 2000), HEAT (Barwick et al. 1997), AMS (Aguilar et al. 2007; Alcaraz et al. 2000) and MASS 
(Grimani et al. 2002). 



siderably improved their model, is still widely used as a 
reference for predicting the positron fraction: we note that 
this model clearly overshoots the AMS data above a few 
GeV, and would therefore underestimate the positron frac- 
tion, should the electron data of AMS be close to the true 
flux and unaffected by systematics or unknown transient 
effects. 

In Fig. 12, we show the positron fraction obtained for 
both the soft (left panel) and hard (right panel, respec- 
tively) electron spectra. For the positrons, we have used 
the Kamae nuclear cross-sections and the Shikaze proton 
and alpha injection spectra. The yellow band is bounded 
from below (above) by the MAX — short dashed curve — 
(MIN — solid curve — , respectively) set of propagation 
parameters, while the central long-dashed curve represents 
the MED configuration. A solar modulation with = 600 
MV has been applied to the positron flux, which corre- 
sponds to the level of solar activity during the data taking 
of AMS. In the same figure, we also report the positron 
fraction obtained with the positron flux of Moskalenko & 
Strong (1998), but with our parameterization of the elec- 
tron flux. 

We see that, in the hard index sizeable excess is 

present in the high energy tail. The MED reference curve 
is marginally compatible with the HEAT and AMS data 
above 10-20 GeV, which instead lies closer to the upper 
border of our predictions, thus favoring the MIN model, 
which is consistent with our predictions of the secondary 



positron flux (cf. Fig. 5), should the positron flux be domi- 
nated by a single secondary contribution. Therefore, when 
the theoretical uncertainties are considered, a clear assess- 
ment of an excess is not statistically significant on the basis 
of the HEAT and AMS data alone, apart from the 2a ten- 
sion with the AMS data point at 12 GeV. Nevertheless, 
in the case of the PAMELA data, the MED reference flux 
is clearly incompatible with the experimental determina- 
tions for energies above 10 GeV. Even when theoretical 
uncertainties in the positron flux are taken into account, 
an excess is probably present for a hard electron spectrum. 

When using the soft electron parameterization instead, 
we see that although an excess might still be apparent, its 
amplitude has strongly decreased, making it of least statis- 
tical relevance. The MED model indeed reproduces all the 
data-sets well from a few GeV up to 40 GeV and a deviation 
is present for the last two bins of PAMELA, where the er- 
ror bars are large due to reduced statistics. The PAMELA 
data may therefore be indicative of an excess also for a soft 
electron spectrum and energies above 50 GeV, but once the 
theoretical uncertainties on secondary positrons and statis- 
tical fluctuations in the data are taken into account, the 
amplitude of this excess is of least relevance. This implies 
that, for a soft electron spectrum, the secondary positron 
yield might still represent a very important contribution 
to the entire cosmic positron flux. We note, however, that 
below 4 GeV, the MED configuration appears to disagree 
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TOA Positron energy [GeV] TOA Positron energy [GeV] 



Fig. 12. Positron fraction as a function of the positron energy, for a soft (left panel) and hard (right panel) electron 
spectrum. Data are taken from CAPRICE (Boezio et al. 2000), HEAT (Barwick et al. 1997), AMS (Aguilar et al. 2007; 
Alcaraz et al. 2000), MASS (Grimani et al. 2002) and PAMELA (Adriani et al. 2008). 



with the HEAT and AMS data, which would favor the MIN 
configuration. 

Thus, we have attempted to discuss the positron frac- 
tion data by considering two different parameterizations for 
the electron flux, both consistent with the data below 100 
GeV. If we had considered other fit parameters that were 
also consistent with the data, by modifying for instance the 
normalizations so as to ensure that they remained corre- 
lated with the spectral indices^ , we would have found quite 
different results opening up more or less the banana shape 
characterizing our predictions in Fig. 12. We therefore em- 
phasize strongly that it is difficult to interpret the origin 
of the positrons observed in the positron fraction without 
comparing secondary positron predictions with the positron 
data, or a precise measurement of the electron flux. This 
requires access to both the electron and positron data sep- 
arately, which have not yet been released in the case of 
PAMELA. With this in mind, it may be potentially unsafe 
to infer strong statements about the possible nature of any 
excess, since its amplitude and shape depend strongly on 
the underlying assumptions, which are not well constrained 
at the moment. 

To conclude this section, we recall that in the previous 
sections we have shown that our theoretical predictions of 
the secondary positron flux perfectly agree with the avail- 
able data (PAMELA data for the electron and positron 
components separately are not yet available), especially 



^ A correlation between the normalization and the spectral 
index obviously appears in the case of a single power law fit, but 
this would have absolutely no physical meaning if there actually 
were different spectral contributions — secondary, far primary, 
local primary — at different energies to the electron flux. A multi- 
component fit would break this correlation. 



when the theoretical uncertainties are properly taken into 
consideration. We have shown in this section that the elec- 
trons play a very important role in the interpretation of the 
positron fraction released by the PAMELA collaboration. 
The scatter in the experimental measurements of the elec- 
tron flux does not allow us, at the moment, to characterize 
fully the amplitude and shape of any positron excess with 
respect to the secondary yield. We have provided a critical 
illustration of this by using two different parameterizations 
of the electron flux {hard and soft cases) that are consistent 
with the data below 100 GeV, from which we have shown 
that the amplitude of an excess might be much smaller than 
expected from the positron/electron model of Moskalenko 
& Strong (1998). Even if small, this excess could originate 
from additional astrophysical processes or dark matter an- 
nihilation. This implies that it is essential to estimate and 
understand the theoretical uncertainties affecting the back- 
ground, as we have attempted in this paper. 

6. Conclusions 

Our aim has been to compute a correct estimation of the 
astrophysical positron flux and the corresponding uncer- 
tainties. 

First, we have compared the various models available for 
the interstellar secondary positron production. It has been 
shown that more positrons are expected when the proton 
flux from Shikaze et al. (2007) is used, as compared to the 
case proposed by Donato et al. (2001) . Moreover, for a given 
proton flux, the three positron production cross-sections we 
have considered produce different results: below a few GeV, 
the parameterization of Badhwar et al. (1977) gives more 
positrons, whereas above a few GeV, the model of Tan & 
Ng (1983) predicts a higher positron production. At any en- 
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crgy, the parameterization of Kamae et al. (2006) produces 
the lowest amount of positrons. 

Concerning the propagation of the positrons in the in- 
terstellar medium, we have used a Green function approach 
that led us to disregard convection and diffusive reacceler- 
ation: we have specifically included diffusion and energy 
losses due to inverse Compton scattering on cosmic mi- 
crowave background photons and synchrotron radiation. 
Nevertheless, this analytical method allowed us to scan our 
~ 1, 600 sets of propagation parameters compatible with 
the boron to carbon ratio measurements, and therefore to 
determine astrophysical uncertainties in the positron flux 
predictions. We showed that varying the diffusion parame- 
ters does not have the same effect as for primary positrons 
(Delahaye et al. 2008). For exotic positrons created in the 
Dark Matter halo, the thickness of the slab 2L was the 
most relevant parameter because the increase of the dif- 
fusion zone implies the increase of the number of sources, 
whereas for sec-ondary positrons - which are created in the 
Galactic disk only - the most relevant parameter is the 
diffusion constant Kq. Therefore, we expect the sets of pa- 
rameters that basically maximize the primary positron flux 
to minimize the flux of secondary positrons, and vice versa. 

We also showed that, because of energy losses during 
propagation, most of the positrons detected at the Earth 
have been created in the nearby 2 kpc: this is the reason 
why we could safely neglect the variation in proton flux in 
the Galaxy. By solving the complete equation Eq. (22) with 
a numerical technique, we proved that all the other effects 
(convection, reacceleration, and other losses) can safely be 
neglected below 10 GeV and that our method is valid. 

Finally, and this is our most important result, our esti- 
mation of the positron flux is compatible with all available 
data. This does not mean that there is no exotic positron 
contribution, since we have not tried to flt the data with a 
single diffusion model. However, this shows that one should 
be cautious before claiming that there is any excess in 
present data. Regarding a possible excess in the positron 
fraction, we have also clearly shown that the electron flux 
plays a role that is as important as that of positrons. This 
might sound tautological because the positron fraction is no 
more than a ratio, but so much energy is involved in sup- 
port to this that we repeat this point. Figure 12 provides 
compelling visual argument. 

The released PAMELA data (Adriani et al. 2008) show 
a clear increase in the positron fraction for energies above 
10 GeV. From our analysis, whose objective is an accurate 
determination of the positron flux, we derive the conclu- 
sion that an excess is clear for a hard electron spectrum, 
while for a soft electron spectrum the rise in the positron 
fraction may be explained by the standard secondary pro- 
duction. By considering in turn the various parameters, we 
find in general that the PAMELA measurements are in 
excess of what a pure secondary component would yield. 
Nevertheless, if the electron spectrum is soft, most of the 
PAMELA data points are aligned with our MED predic- 
tion. We note also that, in that case, the two last energy 
bins feature an increase, but the experimental uncertainties 
are large there and a presence of an excess is, in this case, 
currently not statisticlly significant. 

More insight into these issues will therefore require, 
from the theoretical side, a revised understanding also of 
the electron flux, including the determination of its un- 
certainties, and from the experimental side, the separate 



provision of the elec;tron and positron fluxes, to allow more 
robust comparison of theoretical predictions with the data. 
In addition, the upcoming data on cosmic rays above 10 
GeV will allow us to reduce considerably the theoretical 
uncertainties in all cosmic ray fluxes and help us to eluci- 
date the experimental status of the so-called excess in the 
positron spectrum. 
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